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In nonabelian gauge theories, the Coulomb gauge is one of the most important ones. The theories with this gauge 
are described in terms of physical fields, so that the unitarity is manifest. In the Lagrangian, second-order, formalism 
of Coulomb gauge perturbation theory, there appear, at one-loop order and above, divergent energy integrals of the 
form, 



Here po is the temporal component of the four- momentum P^ L — (po>p) an d indicates a set of external three 

momenta of the amplitude under consideration. [We use capital letter Q M to denote a four vector — (q° , q) and 
q J for denoting a three vector.] 

The following results are established by now on the energy divergence problem. 

1) In the calculation using the Hamiltonian, phase-space, first order, form of Feynman rules, energy-divergences like 

do not appear [jj in the first place. 

2) On the basis of a correspondence formula which equates amplitudes in a covariant gauge to those in a gauge 

without ghosts, Cheng and Tsai [2( showed that, when all relevant contributions are added, cancellation should 
occur between the divergences like (ll.lf) in the Coulomb gauge. 

3) With the aid of an interpolating gauge, which interpolates between a covariant gauge and the Coulomb gauge, it 

was shown in the phase-space formalism that, when all participating contributions are added, cancellation 
occurs between the divergences like (jnj. (See, also, 0.) 

Furthermore, following two types of ill-defined integrals appear: 



The type l|1.2l) appears at one-loop order and above, and the type (|1.3|) does at two-loop order and above. 

In Coulomb-gauge QCD, one encounters a problem of operator ordering in the Hamiltonian. This problem was 
resolved by Christ and Lee 0, along the line of the earlier work of Schwinger Q: The quantum Hamiltonian is 
different from the classical Hamiltonian by special terms, labeled V\ + V?.. Since then, it has been shown that the 
ill-defined integrals of the form l|1.3|l are connected 00 with these V\ + V<i terms, and integrals of the form (|1.2(l can 
be set equal to zero 

As mentioned above, when all relevant energy-diverging contributions are added, energy divergences cancel out, 
provided that the remaining three- momenta integrations are convergent. It has been pointed out [7j that, in dealing 
with the renormalization parts, difficulties arise in simultaneously handling the diverging energy integrations and 
renormalizing the ultraviolet (UV) divergences. A formal proof of cancellation of energy-divergences and algebraic 
renormalizability is given in |3J with the aid of an interpolating gauge in the phase-space formalism. Recently, this 
issue is studied [8j in an example in which quark-loop subgraphs are inserted into the second-order gluon self-energy 



I. INTRODUCTION 
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graphs. As mentioned in 1) above, in the phase-space formalism, two integrals over the internal energies converge, 
provided that the internal three spatial momenta are held fixed. It is found in [8j that, when one first computes each 
subgraphs and performs renormalization, energy-divergences re-appear in the final energy integrals. Thanks to the 
Ward identity, these energy-divergent contributions are cancelled out when all relevant contributions are added. 

In this paper, generalizing the analysis in [8j to all orders of perturbation theory, we give a proof of cancellation 
of energy-divergences without spoiling the renormalizability of the Coulomb-gauge QCD by using the Feynman rules 
derived from the Lagrangian, second-order, formalism. In dealing with the energy-divergences of the form 
extra Vl + V 2 terms do not participate @, and then we can use the Coulomb-gauge Feynman rules derived from the 
Lagrangian formalism. 

In Sec. II, we derive a power-counting formula for divergent energy integrals. In Sec. Ill, we show that the energy 
divergences, Eq. are cancelled out to all orders of perturbation theory. In Sec. IV, for completeness, with the 

help of the power-counting formula obtained in Sec. II, we identify the diagrams that yields ill-defined energy integrals 
(|1.2|) and (|1.3|) . In Sec. V, we show that the renormalization does not spoil the cancellation of energy divergences. 

II. PRELIMINARY 

QCD in the Coulomb gauge is defined, with standard notation, by the effective Lagrangian density 

£ cff = -\fZ v f« v - l-id^id^) 

+d l f) a (S ab d l - gfacbAl) i] b 

+yL (i$ - m - gt a A a) 1> , (2.1) 

where t a = X a /2 and F£ u = d^A v a - d h ' A» - gf abc A^ b l A u c . Here A£ denotes the gluon field and g a (fj a ) denotes the 
(anti)FP-ghost field. We have introduced one quark flavor (ip,ijj). Generalization to the case with several quark 
flavors is straightforward. Generalization to other nonabelian gauge theories is also straightforward. 

Throughout in the sequel, we restrict ourselves to the strict Coulomb gauge (a — > 0). The propagators in the 
Lagrangian formalism may be extracted from the bilinear terms (in fields) of C e s in Eq. 1)2. If) : 

(TAi(x)Al(y)) m iS ab p fj P ^ + =iS ab D»(P), 

(TA° a (x)A° b (y)) i5 ab ^^i5 ab D™{p), 

p 2 

FT 1 ~ 

(Tr] a (x)fj b (y)) -^-4 -i5 ab ^ = i5 ab D(p) . (2.2) 
Here, P 2 = p\ — p 2 (p = |p|), <^(p) = 8^ — p l p 3 jp 2 : , and 'F.T.' stands for taking Fourier transformation. 

Superficial degree of energy divergence 

It is sufficient to deal with one particle irreducible diagrams. ^From them, we take a particular diagram G. We 
introduce the following abbreviations; 

i) S A' for the spatial component of the gluon field, A l a , which we call the transverse gluon (tgluon) in the sequel, 

ii) '0' for the temporal component of the gluon field, ^4", which we call the 'Coulomb', 

iii) 'G' ('&) for the (anti)FP-ghost, 

iv) 'q' ('(f) for the (anti)quark. 

We adopt the following notation to describe the diagram G: 

Ii = number of internal lines of i (i — A, 0, G, q) , 
Ei = number of external lines of i 

(i = A, 0, G, q and q) , 

V$A, V4A, VaAQ, VAQO , VaAOO , VaGG' VAqq, Vogq 

= number of vertices indicated by the suffices . 
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Let us find the superficial degree of energy divergence, ui(G), of G. When u>(G) < —1, energy integral converges. 

To w(G), each loop contributes +1, each internal tgluon line contributes —2, each internal quark line contributes 
— 1, and each l AA0' vertex contributes +1. Then, it is clearly 

w(G) = L — 21 a — I q + V AA0 . (2.3) 

^From the topological structure of G, we have 



L = I a 



2I A + E A = 3V 3A + W 4A + 2V AA0 + V AQ0 

+2VAA00 + V A GG + VAqq , 

2I + Eq = Vaao + 2V A oo + 2V AA oo + Voqq , 

2I q + Eg and q = ZVjiqq + 2Voqq j 

2/ G + £ G = 2V AGG . (2.4) 

In the first equation, summation is taken over all types of vertices in G, where we have used the fact that there is 
energy conservation at each vertex but there is also one overall energy conservation. Using the relations l|2.4|l in Eq. 
(|2.3() . we obtain 

w(G) = 1 + - (E A — Vaao — Vaoo — 2Vaaoo — Vagg) 

-\ (5V 3A + 6V iA + W A qq + V m ) . (2.5) 

The directions of the momenta of the external lines are taken toward outside of the diagram G (c.f., Fig. 1). Eq. 
(|2~5|) tells us that u(G) < 1. 



III. CANCELLATION OF u(G) = 1 ENERGY DIVERGENCES 

A. Participating diagrams 

The diagram G with cu(G) = 1 yields the divergent integral of the form From Eq. I|2.5|l we learn that the 

u>(G) = 1 energy-divergence arises only when the following three conditions are simultaneously met: 

(CI) E A = Vaao + Vaoo + 2Vaaoo + V a q g , 
(C2) E = E G = 0, 

(C3) V 3A = V iA = V A qq = V q q = 0. 

These conditions leads to the following propositions: 
(PI) Energy divergence arises only from the tgluon amplitudes [(C2) and (C3)]. 

(P2) In an energy-divergent tgluon amplitude, the vertices in (CI) above, Vaao, V A oo, V AA qq, and Vaggi are ^he 
external vertices, i.e., external-tgluon lines go out from them. In particular, when one tgluon goes out from a 
vertex VaaoOi if an Y! m the diagram G, no energy divergence arises [(CI)]. 

(P3) Energy divergent diagrams do not have internal vertex [(P2) and (C3)]. 

Then, 

(P4) Then, energy divergence can arise only from tgluon one-loop diagrams, Fig. 1. 
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Fig. 1: One- loop N-tgluon diagram. 



B. Structure of the building blocks of the tgluon one-loop diagrams 



Let us compute the amplitude for the diagram in Fig. 2(a) that is a part of Fig. 1 



A 2a = [-gfbceiP + q° 2 )] 



P 2 + ifH 



^9 fbcefdae 



P 2 + i0 



[gfdae^-p )} 



The second term in the square brackets on the last line does not yield the lo(G) = 1 energy diverg 
and throughout in the following we are concerned only with the portions that participate in the <jj(G) 
divergence. Then, we have 



A 2a 



p l p 3 



ig fbcefdae 5 V — , 



P 



(3.1) 

ence. Here 
= 1 energy 



(3.2) 



where '~' indicates that the right-hand side is the portion of the left-hand side that leads to the energy divergent 
contribution to the one-loop amplitudes under consideration. 



b 



P, e 
(a) 



a 



b a 
(b) 




Fig. 2: Lowest order four-point diagrams, each of which is a part of Fig.l. Here Solid lines represent tgluons; dashed lines 
represent Coulombs; and dot-dashed line represents FP-ghost. 



The amplitude for the diagram in Fig. 2(b) reads 

Alb = —ig 2 {fbcefdae — fbdeface) 
*2b t J^ 2 b 



A 2 b A 2 b ■ (3-3) 



We see from Eqs. (|3.2() and l|3.3[) that the partial cancellation occurs between A 2a and AQ . Similar partial cancellation 

( 2) 

occurs between the contribution from the diagram that is obtained from Fig. 2(a) by (a, i, Qi) <-> (b,j, Q 2 ) and A 2b . 
Thus, we have 

A 2a + 41' - fbcefdae^ (= A 2 ) . (3.4) 
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With understanding that Figs. 2(a) and 2(b) are always combined into the form i|3.4|) . we will forget Fig. 2(b) or 
Vaaoo hereafter. 

The contribution from Fig. 2(c) reads 



A 2 



gfbce(2p J 

^9 fbcefdae o~~ 



[gfdae(2 P i +q[)] 



(3.5) 



Here we have used the fact that, in the strict Coulomb gauge adopted here, q\e 1 r (qi) = (r = 1,2) and q\D lk (Qi) = 0, 
where e^,(qi) and D ik (Q 1 ) are, in respective order, the tgluon polarization vector and the propagator, Eq. 1|2.2|) . which 
are to be attached to A 2c - Thus, q\ may be dropped. Similarly, q\ may be dropped. 
In a similar manner, we have, for the contribution from Fig. 2(d), 



A 2d = 



-gfaceirf - qi) -o [-gfdaeP 1 ] 

J p J 



ig fbce fdae 



p l p 3 



(3.6) 



We observe that, besides the difference between the overall factors, the functional forms of A 2 , A 2c , and A 2 d are the 
same. We note that the integrand of a potentially energy divergent tgluon one-loop amplitude is independent of the 
temporal component of the loop momentum. 

Finally in this subsection, it should be emphasized that the one-loop diagram that includes two or more adjacent 
tgluon propagators does not yield energy divergence (c.f., Eq. jSJJl ), 



C. Absence of overlapping energy divergences 



Here, we show that no overlapping energy divergences arise. Assume that, in the two-loop diagram depicted in Fig. 
3, both the left-side one-loop (S^) and the right-side one-loop (Sr) yield energy divergence. Then, from (P4) above, 
four lines with momenta P, Q, P', and Q' are the tgluon lines. Furthermore, from the observation at the end of Sec. 
IIIB, the line with momentum P — Q is the Aq or Coulomb line. Then, the vertex factor for the vertex at which three 
lines with momenta P, Q, and P — Q meet is proportional to po + <Zo- As can be seen from Eq. (|3.1|l . the 'po' Oqo') 
part participates in the energy divergence for (S^) but does not participate in the energy divergence for S# (5^,). 
Thus, the energy divergence does not arise simultaneously from and S^j. 




Fig. 3: Two-loop tgluon diagram. 



D. Cancellation of u)(G) = 1 energy divergences 

One can read off from Eqs. I|3.1|l - l|3.5() the relative factors for the Coulomb- and the tgluon-propagators; 1 : — 1, 
and the relative vertex factors for a AQQ- and AAO-vertices; 1 : —1/2. To represent these relative factors, we introduce 
matrices, 
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Here the first rows and columns correspond to l A' and the second rows and columns to '0'. Vn = comes from the 

(G) 

fact remarked at the end of Sec. 111B. Let A N be the sum of TV-point tgluon amplitudes for Fig. 1, where the loop 
consists of Coulomb and/or tgluons lines. Then we obtain for the w(G) 



(G) 

1 energy-divergent contribution to A N , 



A 



(G) 

N 



Tr 



JV 



A 



(GO) 
JV 



(3.8) 



where A^°^ is the contribution from Fig. 1, where all N propagators are the Coulomb ones. Note that, for N 
symmetry factor 1/2 is necessary, which is included in A N and A N . Through mathematical induction, we obtain 



N 

VV) = 2 



-N 



so that 



A 



(G) 
N 



-N+l N 
—N N+l 



)1-A^(G0) 



JV 



(3.9) 



(3.10) 



The tgluon amplitude A N for the FP-ghost one-loop diagrams is obtained from through the following 

operations (cf. Eqs. I|3.5|l and J^JJ); a) change the sign of each propagator, b) multiply a factor —1/2 for each vertex, 

c) multiply a factor 2, for 3 < N, that corresponds to two diagrams with opposite circulation of fermion number, and 

d) multiply —1 that comes from one fermion loop. Thus, we obtain 



A 



(FP) 
N 



R-2-R 



TV 



■A 



(GO) 
N 



yl-N A (G0) 



N 



(GO) 



For 3 < N, the factor '2' in the mid-term comes from c) above. For N — 2, as mentioned after Eq. (|3.8|l . A 
includes a symmetry factor '1/2', while, A^ P ^ does not include a factor '1/2'. Then, the factor '2' is necessary in the 

(G) ( P PI 

mid-term. Thus, we see that A N cancels ^4 ' 



A 



N 



A^j ~\~ A 



(FP) 
N 



(3.11) 



E. Regularization 

As mentioned above, A^ ^ is independent of the temporal component of the loop momentum and is of the form 

A<£ 0) = Jd 3 P J dp F(p,...). (3.12) 

Integration over p$ diverges. In order to rigorously handle this integral we should introduce some regularization. It 
is desirable to adopt the regularization that preserves the BRST symmetry of the theory, so that the regularizing 
quantum effective action T satisfies the Zinn- Justin equation, whose explicit form is not necessary for our purpose. 
To our best knowledge, two candidates are available, i.e., the interpolating gauges 0, and the split dimensional 
regularization 

Throughout in this paper, we adopt the interpolating gauge 3], which is defined by introducing the gauge condition, 

£3,41 + Mi = o 

(£ is a real parameter). The Coulomb gauge is obtained by taking the limit £ — » 0. £ here plays a role of regularization 
of the energy-divergent integrals in the Coulomb gauge. For arbitrary £, T satisfies the Zinn- Justin equation, which, 
in the limit £ — > 0, turns out to be its Coulomb-gauge counterpart. 
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Gluon- and ghost-propagators read 

A y '(P) 

A 00 (P) 
A 0l (P) 
A(P) 



1 



P 2 + ^0^ 



P 2 -C(2-Qpg 
(P-P c + i0+) 2 ^ 



P 



(P • P f + ?:o+) 2 ' 

A <0 ^OP* 



(P-Pt + i0+) 2 



1 



P • P 4 + W 



A U0 (P) , 



(3.13) 



where P^ = (£po,p)- 

In the following, we only keep the terms that turn out to the w = 1 energy-divergent ones, and we ignore the terms 
that tend to in the limit £ — > 0. In place of Eqs. (|3.4() . I|3.5(l . and (|3.6|l we have, in respective order, 



A 2 
A 2c 



■ a , , P 2 -£(2-flpg 4 , 

"«5 fbcefdae /ri „ , . n ,-> V 



4i.g fbcefdae 



A 2 d — fbcefdat 



(P • Pf + i0+) 5 

o 

(pTp7+ z0+) 2 



p'p 3 



(P • P e + «0+) 



* 7 

rPF 



(3.14) 
(3.15) 
(3.16) 



For obtaining .A 2c for Fig. 2(c), we have used the fact that gje*(qi) = O(0 (r = 1,2) and <7jL> lfc (Qi) = 0(0, 
where e*(qi) and D lk (Q 1 ) are, in respective order, the tgluon polarization vector and the propagator, which are to be 
attached to Aic- Then, q\ that was present in Eq. I|3.15|l (cf. Eq. I|3.5|l ) turn out to be of 0(0 and lead to vanishing 
contribution in the limit £ — > 0. Similarly, q\ that was present in Eqs. (|3.15(l and l|3.16|l can be ignored. 



b 

Q2V 



a 



P, e 
(a) 



N d 



P,e 
(b) 



Fig. 4: Lowest order four-point diagrams, each of which is a part of Fig.l. 
In the interpolating gauge, two diagrams as depicted in Fig. 4 also participate. Corresponding amplitudes read 

ifPoP* 



Aia = [gfbce(2p - q 2 ) 3 ] 



-2ig fbcefda 



(P-P e + iO+) 2 



[gfdaeiq !-?")} 



(P • P e + ^0^ 



46 • 



Diagrams as depicted in Fig. 5 do not yield energy divergence in the limit £ — ► 0. 



b 



P,e 
(a) 



a 

^ d 



b 



a 



P,e 
(b) 



Fig. 5: Lowest order four-point diagrams. 
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Let us make a change of variable, po — > Po/V£- Then, ignoring the terms that lead to vanishing contribution in the 
limit £ — > 0, we have 



A 2 
A 2c 



-i9 2 fbcefdae(p 2 ~2p~ 2 )V^(P) 
^fbcefdaeP^iP), 

A 2d * ig 2 f bce f dae (p 2 ~p 2 )V l \P), 
A ia =s A i b^-2ig 2 f bce f d acf 'D ii {P), 



(3.17) 
(3.18) 
(3.19) 
(3.20) 



where P = (po,p) and 



V ij {P) 



p % pj 



(P 2 + i0+) 2 

^From these equations, we obtain, in place of Eq. I|3.7|) . for V and V, 



V l = 



2p 2 Q /p 2 - 1 p 2 Jp 2 



1 



, v = 



-1/2 
-1/2 1 



(3.21) 



The contributions that turn out to ui(G) = 1 energy-divergent contributions in the limit £ — > is given by (cf. Eq. 
(E3), 



4 ( G ) ~ _L f dp° 



where A^ ^ = (2ir^/£) 1 J dpoA^' 1 . Through mathematical induction, one can show that 



.(GO) 



VI J 2^ 



Tr 



' N 

n H 



U=l 



^ G0) + o(e /2 ) 



(3.22) 



Tr 



' AT 

n (w) 

./=i 



iV 



= 2 



N 2 _ ~2 
1-JV J~J Pi ~ PjQ 



i=l 



P/ 



Then, Eq. (|3~2ll turns out to 



A 



(G) 
N 



>1-JV 



v 7 ? 7 2vr 



No ~ 2 



1=1 



Pi 



A^+oie' 2 ). 



(3.23) 



The amplitude for the FP-ghost one-loop diagram is obtained similarly as in Sec. HID using Eq. 13.13fl or Eq. 
(ETUI : 



A 



(FP) 
N 



)1-JV 



dp 



JV 



2 ~2 

Pf ~ Pfo 



n * 



A^+o(e /2 )- 



(3.24) 



Then, we see that 0(£ -1 / 2 ) contributions to Affl and A^ P ^ are cancelled out, and, in the limit £ — > 0, we have 

Thus, Eq. (|3.11(l gets a sound foundation. 

It should be emphasized that the diagrams in Fig. 4, which are absent in the strict Coulomb gauge, participate 
here. 

Analysis with the split dimensional regularization leads to the same result, which we do not reproduce. 



IV. ILL-DEFINED INTEGRALS WITH lu(G) = 

In this section, for completeness, we briefly mention the diagrams with lo(G) = 0. From Eq. I|2.5|) together with 
the observation in the last section, we see that oj(G) = energy divergences and ill-defined integrals arise from the 
following diagrams (see, also, Q): 



Fig. 6: One-loop g-g-tgluon(s) diagram. The double line 
represents a quark. Here each solid line that constitutes 
the loop is a Coulomb or a tgluon. 



Fig. 7: Two-loop tgluon diagram. The double line repre- 
sents a quark, each solid line that constitutes the loop is 
a Coulomb or a tgluon. 



1) Gluon one-loop amplitudes, of which a number of external Coulomb field is at most one. 

2) g-g-tgluon one-loop amplitudes of the type as shown in Fig. 6. 

3) Tgluon two- loop amplitudes (Fig. 3). 

4) Two-loop amplitudes of the type as depicted in Fig. 7. 

Now let us inspect Fig. 3. Cutting the line with momentum P (or Q or P — Q), we obtain an one-loop diagram. If 
u> = 1 energy divergence arises from this one-loop diagram, as has been proved in the previous section, it is cancelled 
out together with its relatives. Then, from Fig. 3, integrals of the forms 1|1.2|) and (|1.3|l emerge, which are ill-defined 
ones. The integrals of the form l|1.2f) emerge from 1) - 3), and the integrals of the form (|1.3|) emerge from 3) and 4). 
As mentioned in Sec. I, the form (II. 2|) can be set equal to zero, and the integrals of the form (|1.3|l are connected with 
V\ + V2 terms of Christ and Lee |4| . 

This is, however, not the end of the story. As a matter of fact, newcomers enter into the stage through renormal- 
ization, i.e., the renormalization counterterms. In the next section, we deal with this issue. 



V. COMPATIBILITY OF CANCELLATION OF ENERGY DIVERGENCES AND 

RENORMALIZ ABILITY 

We adopt the minimal subtraction scheme in dimensional regularization and introduce the counter Lagrangian 
density 5C c g that includes different renormalization constants. The effective QCD Lagrangian density £ c ff with 
addition of 5£ e g in the strict Coulomb gauge reads 01 

£ cff + 5C cS = -^f {d^A^Ai - ihM,;),X) + Z 32 d a A a d z A° a - 

-^-(diA^Ai) 
la 

+gZ ll f abc (d l A a J )A\Al + gZ 31 f abc d Q A l a AlA a + gZ 32 f abc d l A a A H A° c 
2 2 

— ^-Z^fabcfadeAbiAcjAjAi — Z 3 ij ' abc j ' a deA%A ci A d A\ 

+Z 3 (dif] a )d % r] a ~ gfabcidtfj^Alrjc 

+Z 2 tp {1$ - (m - dm)) i> - gZi&aAl^ - gZ 2 ^t a A a [ 1 ^ , (5.1) 

where the limit a is understood to be taken. Z's and Z's in Eq. (15. 1|) obey the Slavnov- Taylor identities in the 
narrow sense, 

Z 3 i _ £ _ Z\\ _ Z 32 _ Z 33 _ Z 2 (5 2) 

Z\\ Zn Z 3 i Z 32 Z\ 

The counter Lagrangian 5£ c ff introduces new vertices, which we call counter vertices. Repeating the argument in Sees. 
II and III on the superficial degree of energy divergence, we find that the energy divergence arises from the diagrams 
that are obtained from each energy-diverging diagram Gi (i — 1,2,...) treated so far by inserting the counter vertices 
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in all possible ways. Let A[^ (j — 1,2,...) be the amplitude for thus obtained jth diagram, which may be computed 
from the Lagrangian density C e s + <5£ c ff in Eq. (|5.1|l . Summing over all Wf^'s, we obtain A — J2i + J2j -^J i 

where A\ is the amplitude for Gi. 

We now show that A is free from energy divergence. Eq. (|5.1|) tells us that the propagators and the vertex factors 
extracted from C c a + SC c a arc obtained from their bare counterparts through the following replacements: 



D 00 (p) — * Z 33 1 D 00 (p) , 





♦ Z 3 l D{p) , 


yijk _ 
abc 


7 \rijk 


yijO _ 
abc 


"* ^31 Kibe ' 


yiOO _ 
v abc 


' abc i 


yijkl 
abed 




yijQQ _ 

abed 


Z 31 ^afced 


v\ - 

abc 


-> v\ • 

abc 







(5.3) 

where V£ bc is the three-tgluon vertex factor, V % J bc is the A l a A 3 b A^-veitex factor, V* bc is the fj a A l b -q c - vertex factor, 
and so on. Explicit form of the quark propagator is not necessary for our purpose. 

We first show that the key equations (|3.4|l and i|3.11[l remain to hold. From Eq. I|5.3[) . we see that Aia-, Eq. (|3.2[) . 
and A 2b , Eq. , turn out to be 



A 2a 



Z 2 

~7i — A 2a = Z-sxAla , 
^31 



and 



A 



ZsiA^ b , 



respectively. Then, the same partial cancellation as for Aia and A^ b in Sec. Ill occurs, and we obtain, in place of 
Eq. (TOl . 



A 2a + A {1) 



21, 



Z31 ^A.2a 



A 



Z31A2 



(5.4) 



Taking A^ ^ as the reference amplitude, we have, in place of the matrices V and V in Eq. l|8.7[l . 



V = 

V = 



-Z^Zm 
1 





-Z 3 iZ 32 /2 

Using the identities (|5.2p . we have, in place of Eq. Q3.9[l . 



-■^31-^32 /2 

1 



iV 



^PVJ = 2 

Then, we see that the relation (|3.10(l holds as it is 

A 



-N + 1 N(Z 32 /Z 31 ) 
-N(Z 31 /Z 32 ) N + l 



(G) 
N 



N 



The tgluon amplitude for the FP-ghost one-loop diagrams is obtained using Eq. (|5.3|) . 



A 



(FP) 

■N 



,1-N 



Z 3 



■^33 
Z32 



N 



A 



(GO) 



N 



where use has been made of Eq. (|5.2() . Thus, A N cancels A N 



A 



(G) 
N 



A 



(FP) 
N 



(5.5) 
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Proof of Eq. 15.51) with the aid of the regularization by an interpolating gauge 

The effective QCD Lagrangian density with addition of the counter Lagrangian density in the interpolating gauge, 
£,d A° a + diA\ = 0, reads 



£ e ff + SCcff 



Z 3 1 
1 

2^ 



(d„A aj d»Ai - n,M,n,.\:,) + z 32 d Aid i A° a - ^d t A a d l A a 



+gZ 11 f abc (d l A a J )A\A{ + gZ 12 f abc d°A a A° b A ci + gZ 13 f abc d*A a A bl A° c 

g 2 <r 

— -^Z4if abc f a d e A bi A c jA Z d Al — '—Z±2fabcfadeA b A ci A d A l e 

+Z 31 (d l fj a )d l i] a + £Z 32 (dof]a)d r) a - gf a hc{d ll 'n a )A^'q c 

+Z 2 ip (i$ - (to - 5m)) ip - gZ^n^taA 3 ^^ - gZ. 4 , 12 ipt a A°<y '4> , 

where = (£<9q, V) and a — ► is understood to be taken. Z's and Z's obey 



Z31 

~z7i 

Z31 

Z\2 



— Z32 — 



Zja 


Z32 


Z33 


Za 


Z\2 


~ z^ ' 


Z\2 


Z32 


z 2 


Z&2 


Z\3 


Z^\2 



Z-2 



Z^w 
(= D') . 



(5.6) 

(5.7) 
(5.8) 



These equations are derived in a similar manner as in 10]. A few comments are in order. 

• Diagrammatic analysis shows that fjA^rj three-point function is UV finite, as it does in the Landau gauge. 

• In the strict Coulomb-gauge limit, £ — ► 0, D' = 1. 

Gluon- and ghost- propagators can be read off from Eq. I|5.6|l : 



1 



1 



A 00 

A 0i = A 
A = 



z 31 p 2 + ?:o H 
1 1 
z 31 p 2 + W 
Z31 



8 lJ 
~P 2 



p 2 - 2e(z 32 /z 33 )pg + e(Z3i/z 3 3)p 2 o 

(p 2 -(Z 31 /Z 32 )frt+iO+) 2 

p 2 -2tj(D>/D)p 2 +e(D , /D) 2 p 2 o , i 
(p 2 -(D'/D)& 2 + iO+) 2 PP 



p l p> 



Z 2 2(p 2 -(Di/D)Zp 2 +i0+) 2 

i0 _ Z 31 &0P l 



Z 2 2(p 2 -(D>/D)S;p 2 + i0+) 2 



1 



1 



1 



iZ 32 p 2 - Z 3lP 2 + ^0+ Z 3X £(Z 32 /Z 31 ) P 2 ~p 2 + i0+ 

1 {D'/D)t,pl-p 2 
Z 3l (P 2 - (D'/D)ip 2 - z0+) 2 ' 

where use has been made of Eqs. (|5.7|l and l|5.8|l . Through the same manner as above leading to Eq. 
obtain for V; and V, 



(5.9) 
iPHTJ) . we 



Vi 
V 



(D/D') 2 [2(D'/D)p 2 /p 2 -l] p 2 Jp 2 
Pl /Pi I 
-Z 12 /{2Z 13 ] 
-Z 12 /{2Z 13 ) 1 

-D'/{2D) 
-D'/(2D) 1 



where pio = y/^pio- Matrix multiplication yields 



p ^_l(A i -l D/D> 
' 1 2 \ -D'/D A, + 1 
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where 



Mathematical induction yields 



JY 



v (±) 



Then, in place of Eq. (|3.23() , we obtain 



A, = 1 - (D'/D)p%/p 2 



n (a#) = ± 



a 



N 



(+) 
IN <% 

JV 



.4 



(G) 



jl-iV 



2tt 



7iv 



D 7 



i=l 



W 1 

si, 



n 



p1-{D'/D)Pl 

pf 



A 



(GO) 
N 



Using Eqs. (|5.6[) - (|5.8|) and (|5.9|) . we obtain, in place of Eq. I|3.24[) . 



A 



FP) 
N 



yl-N 



dpo 1 
V? J ^ zg 

Z13 



pf - (D'/D)p. 



' N 

n 

.1=1 

N 

A (GO) 
N 
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pi 



il-N 



dpo 



n 

U=l 



pf-(D'/D)p, 



Pi 



A 



(GO) 
N 



Then cancellation occurs between ^4^ and A N , and then, removing the regulator £ — > 0, we have 



(G) , AFP) 



TV 



Thus, Eq. (|5.5(l gets a sound foundation. 



Compatibility of cancellation of energy divergences and renormalizability 

What we have shown above is that, when the renormalization counterterms are included to all orders of perturbation 
theory, energy-divergent contributions are cancelled out. Then, by expanding Z's and Z's in powers of g 2 , energy- 
divergent contributions are cancelled order by order in perturbation theory. 

Armed with this proposition, we are now in a position to show that the cancellation of energy divergences are 
compatible with the renormalizability of the theory. 

We are concerned with the tgluon vV-point one- loop amplitudes A N (N = 2,3,...). Let be the 0{g N+M ) 

contribution to An and Gi (i = 1,2,...) be the set of one- loop diagrams that contributes to A^ ■ For regularizing the 
energy divergence, we employ the interpolating gauge, and to handle the UV divergence we employ the dimensional 
regularization by continuing the spacetime dimensionality from 4 to d. Then, for £ 7^ and d 5^ 4, all contributions 
to An are free from energy- and UV-divergences. 

(I) 0(g N ): As shown in Sec. Ill, A$ is free from energy divergence, i.e., finite in the limit £ — > 0. For 5 < N, A^ 

is finite in the limit d — > 4 (UV finite). For 2 < N < 4, A^P is written in the form, 

Algebraic renormalizability (for arbitrary £) described above in conjunction with Eqs. (|5.6|) - l|5.8|) indicates that 

and are finite in the limit d — > 4 (UV finite). Furthermore, thanks to the above proof of cancellation of 

energy divergences, we have J-"jS = 0, so that, by removing the regulator £ — > 0, we see that Ajj (= G^) is f ree 
from energy divergence and UV finite. 
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(II) 0(g N+2 )\ Through inserting a single one-loop renormalization part into each diagram Gi in all possible ways, 
we obtain a set of diagrams G^ (j — 1, 2, ...). Let (a^^ be the amplitude for G-^: A$ = J2i j (^tv^) 



As mentioned in Sec. IIIC, no overlapping energy divergence arises in ( A N j . Energy divergences and/or UV 

(2) 

divergences that arise from one-loop subdiagrams have already been managed at the first stage (I). Then, A N 
may be written in the form, 



^ = ^ 2) + ^ 2) + o(^ 2) 



The same argument as above in (I) leads .4^ Qpp to be UV finite. 
(Ill) Higher orders: With the aid of the Bogoliubov-Parasiuk-Hepp-Zimmerman (BPHZ) prescription ^lj, one can 
proceed to higher orders and verify that A^ 1 ^ is free from energy divergence and UV finite. 
Finally, it is worth mentioning that, in the course of perturbative computation, if infrared and/or mass singularities 
arise, we introduce small mass for gluons. It is well known that, in any reaction rate, the cancellation occur between 
them El. 
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